Interaction of N solitons in the massive Thirring model and 
optical gap system: the complex Toda chain model 

V. S. ShchesnovichQ 

Department of Mathematics and Applied Mathematics, 
University of Cape Town, Private Bag 7701 Rondebosch, South Africa 

(Dated: February 8, 2008) 

Abstract 

Using the Karpman-Solov'ev quasiparticle approach for sohton-sohton interaction I show that 
the train propagation of well separated solitons of the massive Thirring model is described by 
the complex Toda chain with N nodes. For the optical gap system a generalised (non-integrable) 
complex Toda chain is derived for description of the train propagation of well separated gap solitons. 
These results are in favor of the recently proposed conjecture of universality of the complex Toda 
chain. 
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I. INTRODUCTION 



Recently the complex Toda chain attracted much attention as a possible candidate for 
description of the pulse interactions in integrable and non-integrable nonlinear evolution 
equations [Q, 0, ||, ||, ^ ||, ^ |T0[- For instance, it was shown that the complex 
Toda chain describes the soliton train propagation for all the nonlinear evolution equations 
associated with the NLS hierarchy 0. Quite recently the complex Toda chain was derived 
for the modified NLS equation [|I^, an integrable generalisation of the NLS equation, which 
is associated with the quadratic bundle. 

The complex Toda chain is an integrable generalisation of the well-known real Toda chain 
(see, for instance, Refs. [|, H). In Refs. |P, and [|l^ the comparison of the complex 
Toda chain predictions with the numerical solutions of the NLS and MNLS equations has 
been performed and a good agreement has been established for various choices of the initial 
parameters of the solitons in the train. 

It is noted that the complex Toda chain arises as an approximation of the evolution 
equations describing the inter-pulse interaction in the train comprised of well separated 
solitons with nearly equal amphtudes and velocities. The exponent of the (negative) sepa- 
ration between the solitons serves as the small parameter for the asymptotic expansion and 
derivation of the complex Toda chain can be based either on the variational approach (see 
Ref. 0) or on the adiabatic perturbation theory for solitons (see, for instance, Refs. P, p^). 
However, as noted in Ref. []1T| the variational approach should be used with care. The ap- 
proach based on the adiabatic perturbation theory is equivalent to the Karpman-Solov'ev 
quasiparticle method for the two-soliton interactions [T^. This approach was developed in 
Refs. 

If the nonlinear PDE is not integrable but possesses stable soliton solutions, then the 
train propagation of solitons is described by a generalised (non-integrable) complex Toda 
chain as it is pointed out in Ref. p. 

The complex Toda chain allows a rich class of asymptotic regimes of the soliton train 
propagation @, §], 0]: asymptotically free propagation of solitons, ii) A^-soliton bound 
states with the possibility of a quasi-equidistant propagation, in) mixed asymptotic regimes 
when part of the solitons form bound state(s) and the rest separate from them, iv) regimes 
corresponding to the degenerate and singular solutions of the complex Toda chain. The 
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rich variety of dynamical regimes of the complex Toda chain indicates that it is a good 
candidate for analytical study of the soliton trains. Here I should point out that only few 



simple regimes are exhibited by the real Toda chain |I3| , , thus it is essential to have the 
complex Toda chain in description of the soliton trains. Moreover, the phase space of the 
complex Toda chain with nodes is 4 A^- dimensional, which is precisely the number of the 
real parameters in the train of solitons. 

In the present paper I consider the A^-soliton train propagation governed by two intimately 
related nonlinear PDEs, one of which is integrable and the other is not: the massive Thirring 
model of the classical field theory |]15|, |l^ and the optical gap system Jl^, |18[ [l^, |2^, 



23, 24, 25, 26, 27, 28, 29, 30 



For the massive Thirring model I show that the train propagation of well separated soli- 
tons with nearly equal amplitudes and rapidities is governed by the complex Toda chain. 
Moreover, I derive a non-integrable generalisation of the complex Toda chain, which de- 
scribes the train propagation of well separated gap solitons with nearly equal amplitudes 
and velocities. 

The gap soliton propagation through a grated optical fiber was manifested in recent 
experiments pl|. 



3^ . The results of are of particular interest: there the multiple 



gap soliton formation was observed. 

Localised solutions in nonlinear media with periodic band gaps have a great potential for 
technological applications. One of the most important band gap structures in optics is given 
by an optical fiber with periodic index grating along the axis. From the Floquet-Bloch theory 
of wave propagation in periodic structures it is known that there are forbidden frequency 
bands or band gaps for linear waves. On the other hand, nonlinear wave propagation in 
such structures is possible for the central frequency of the wave packet lying in the band 
gap. Such nonlinear wave is usually called gap soliton. 

The optical gap system was derived within the coupled mode approach for nonlinear wave 
propagation in optical fibers with grating (see, for instance, Ref. [^). It reads 



{EiT — Eix) + i?2 + ( 



0, 



(1) 



i {E2T + E2x) + Ei + + pl^aT) E2 = 0, 

where Ei and E2 are the slowly varying envelopes of two counter propagating waves coupled 



through the Bragg scattering induced by the grating (the hnear cross-couphng terms), the 
nonhnear terms account for the self- and cross-phase modulation effects. The parameter p 
(p > 0) at the self-phase modulation term may range up to infinity [Q, in which case the 
optical gap system models dynamics in the nonlinear dual-core asymmetric coupler p7|. 
Setting p = in the system (|ip one obtains the massive Thirring model. 

The gap solitons in optical fibers were theoretically studied in many works |T3, |TU|, 120, 

1^ (see also the latest review Ref. p9|). Recently, the relevance of 



21|, |2|, ^, ^, ^, ^ 

the system (|^) for description of the optical gap solitons was analytically and numerically 
validated . The general family of the gap solitons was derived in [jl9| using the similarity 
with the massive Thirring model. Recently it was shown that the gap soliton becomes 



unstable when its amplitude grows above some fixed value |]36| , |37 |. 

The soliton-like solutions, which are similar to the gap solitons, were found in nonlinear 
diatomic lattices ^ ^ and in the quadratic (x''^^) materials with a spatially periodic 
linear susceptibihty (grating) |§T],|2|,|3|,|4[|5],|6|,|^|8l. Also, it was shown [|0|, ^ 



that, at some limit, the equations governing nonlinear wave propagation in quadratic media 
with grating and in diatomic lattices are similar to the system (^ though the underlying 
physics is different. 

In the next section, section 2, I state the main results of the paper on the soliton train 
propagation for the massive Thirring model and optical gap system. The details of the 
derivation are placed in the following sections: section 3 for the massive Thirring model and 
section 4 for the optical gap system. The last section contains discussion of the results and 
suggestions for further work. 

II. MAIN RESULTS 

Before formulate the main results I would like to remind some facts about the models 
under study. Let us begin with the massive Thirring model of the classical field theory 
r|,|T§: 



i ivt — Vx) + u + \u\ V = 0, 



(2) 



i {ut + Ux) + V + \v\'^u = 0, 
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where u and v are complex variables, t and x are the time and space co-ordinates, respec- 
tively. The system (Q) is Lorentz-invariant: 

X — tanh(7/)t ^ t — tanh(y)x 



; 1 - tanh^y ) 1/2 ' ( 1 _ tanh2?/) 1/2 ^ 

with M and v transforming as components of the Lorentz spinor. 



In the relativistic kinematics the parameter y is called "rapidity". (Rapidities of two con- 
secutive Lorentz transformations simply add together.) 



The massive Thirring model is integrable by the inverse scattering transform method [|16 
For instance, its one-soliton solution can be written as 

i sin(2'i9) exp{-y/2 + iO) 



u 



cosh(z — i^) ' 

(3) 

zsin(2t9)exp(?//2 + ze) 
cosh(2; + I'd) ' 

where < < 7r/2 and 

z = sin(2^9)cosh(?/)(xo(t) - x), 6 = - cot(2?9)tanh(y)z + (4) 

The soliton has four independent real parameters: y, Xo, and 6. The first two give the 
soliton amplitude and rapidity, while the rest two are the soliton position and central phase 
(the phase at a; = Xo), respectively. The position and phase parameters depend on time: 

dx d(5 

— ^ = tanhy, ~ ~ cos(2'(9)sech?/. (5) 

The first equation defines the soliton velocity: V = tanhy. 

By a suitable Lorentz transformation the rapidity of a Thirring soliton can be put equal 
to zero and the soliton solution (|^) reduces to the quiescent soliton of the massive Thirring 
model. 

The following ansatz is called the soliton train 

JX, i sin(2t9Q,) exp {-yaf^ + i&a) 

V = y 

(6) 

^ _ y> i sin(2^a) exp {ya/2 + iOg) 



where and 0^ are given by formulae similar to (^) and @. It should be stressed that, for 
each soliton, all four soliton parameters in formula (|^) are considered to be t-dependent. 

The CTC for the MTM soliton train. Assume that the A^-soliton train given by 
(I) consists of well separated pulses with nearly equal amplitudes -Qa and rapidities ?/q,, 
numerated by a = 1, in such a way that Xa^\ — Xq, > (here and below denotes the 
position parameter "xq" for the a-th soliton) . Mathematically these conditions are expressed 
as 

|^?a-?9|<?9, |?/„-y|<l, |x„-x„±i|>l, 

I sin(27?o) cosh?/a — sin(2'(9) coshyHxQ, — Xa±\ \ <^ 1- (7) 
Here (and throughout the paper) 'Q and y denote the averages: 

_ 1 TV 1 ^ 

^ = ]7E^- y = j^Y.vo.. (8) 

Define the following new variables: a modified time 

r = sin(2t9)sech(y)t, (9) 

an average phase 

5 = -cos(2??)sech(y)t, (10) 
and the following complex variable for each soliton 

qa = — sin (2'^) cosh (?/)a;Q, — i[5a — 5 — cos(2'(9)sinh(?/)xQ, + avr] + 2a ln[2 sin(2'i9)]. (11) 
Then in the first order of the soliton overlap parameter e, 

e ~ exp{-| sin(2^9^) cosh(?/„)xc, - sin(279«±i) cosh(?/c,±i)a;a±i|}, (12) 
the following two statements are claimed: 

1. the average values d and y do not depend on t; 

2. evolution of the quantities go,, a = 1, ■■■,N, is given by the complex Toda chain with 

nodes: 

^ = _ a=l,...,N, (13) 

dr^ 

where Re{go} = oo and Re{gAf+i} = — oo (i. e., xq = — oo and xn+i = oo, see (pA|)). 



The set of inequalities (0) is similar to the inequalities for the NLS soliton train in Ref. [Q. 
Now I will formulate similar result for the train propagation of pulses governed by the 



optical gap system [2S 



t {Eu - Ei^) +E2+ (1^2 r + pI^iT) E, = 0, 

(14) 

I {E2t + E2,) + Ei + (1^1 12 + p\E2\^) E2 = 0, 



The soliton solution of the optical gap system (14) moving with the velocity V = tanh?/o 

/ 



reads [19 



Ei{x,t) \ e^'^^^'*) 




(15) 



E2{x,t) i [l + pcosh(2?/o)]V2 

where v and u have the form of a Thirring soliton, i.e., given by formulae (@)-® (with 
y ^ Ho)] the additional (nonlinear) phase ip is 

if) = — ^ arctan(tan-i^tanhz), (16) 

1 + pcosh(2?/o) 

with z as in equation (^. 

Here it should be pointed out that the gap soliton becomes unstable when the soliton 
amplitude grows above certain threshold {{^thr ~ '/r/4, see Ref. |36|). This scenario is also 



possible for the train of gap solitons. This instability is the result of the soliton-radiation 
interaction and is beyond the scope of the adiabatic approach. However, being interested in 
stable gap solitons, one can impose the condition i^q, < i^^thr, for all a = 1, ...,N. 

The ansatz I use for the train of well separated gap solitons is given by application of 
the transformation (^) to the train of well separated Thirring solitons (in this case yo = y)- 
Due to the inequalities (|^), the additional phase tpa of each soliton in the train can be 
approximated by formula (|16D with {} = {} and Za = sin(2^9)cosh(?/)(xa — x). 

The generalised CTC for the train of gap solitons. Assume that the train of N gap 
solitons consists of well separated pulses with nearly equal amplitudes and rapidities ya 
numerated by a = 1, in such a way that Xa+i — Xa > and that the conditions (^ are 
satisfied. Associate the following variables with each gap soliton 

Qa = — sin (2'(9) cosh (|/)xq, — i{6a — 6 — [cos(2'i9) — psin(27?)(?/a — l/)]sinh(?/)a;o + an} 

+ 2aln[2sin(2^)], (17) 
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where 



^ 4ptanh(2y) - 
p + sech(2|/) 



Define the modified time r and average phase 5 as in formulae @) and (|lOl). Then, in the 
first order of the sohton overlap parameter e ([l^), the following is claimed: 

1. the average values d and y do not depend on t; 

2. evolution of the quantities Qa, a = 1,...,N, is given by the following generalised 
complex Toda chain with nodes: 

^ = (1 + Ap)(e<3'^+^-«" - e^'^-Q"-) + fi^(e«^+i-«S _ e«^-«^-i), (18) 
where Re{Qo} = oo and RejQAr+i} = — oo. 

Equation (|T8|) is valid for arbitrary values of the self-phase modulation parameter p. 
Here Ap and Bp are p-dependent coefficients: 



with 



Ap = ^{i^ - K,p + i[K{l + u) + p]}, 



Bp = + Kp-i[K{l + u) - p]}, 



ptanh(2y) - sin(4^9) _ 4p(279 cot(2^9) - 1) 



(19) 



p + sech(2y) sin2(2'(9) ' " p + sech(2?/) 

Setting p = in (p!7| ) and (|T8|) one obtains the complex Toda chain for the soliton train 
of the massive Thirring model. 

III. THE CTC FOR THE MTM SOLITON TRAIN 

I will use the adiabatic perturbation theory for derivation of the complex Toda chain. 
Recently the perturbation theory based on the Riemann-Hilbert problem was developed for 



the solitons of the massive Thirring model [49|. For instance, we have derived the following 



evolution equations for the soliton parameters in the adiabatic approximation: 

1 7 dz 



I U9 ^1 777^,^^{e'[r±iki,z) + rl{kl,-z)]}, (20) 

J cosh 2^; + cos 2^ 



dt 2cosh?/ J cosh(22;) + cos(2'i9) 



dy 1 f dz 



dt coshy J cosh(22;) + cos(2'i9) 



lm{e'[r^{h,z)+rl{kl,-z)]-2r\\{z)}, (21) 



^ = tanhy , , (22) 

dt ^ 2sin2(2?9)cosh^?/ ^ ^ 

- = - cosf2^)sechy - M^} + tanhy cot(2^)Re{J} 

dt V ; y 2sin(279)cosh2/ ' ^ ^ 



Here 



d 

J = / — — — - 2|2;z[cos(2^9) + i sm(2^) tanh y] + 

J cosh(2z) + 003(2^) ^ 

— oo 

+ {e^[cos(279)(l - 2z) - 2izsm{2^) tanhy] + e-'}[r^{ki, z) - rl(A;*, -z)]). (24) 
The functions r± and ry are given by 

r,(k,,z,t)^ — [h--K^-y (25) 

where /ci = — exp{— y/2 — id}. The "variational" derivatives denote fictitious evolution of 
the u and v as if under the action of the perturbation only (in other words, either i6u/6t or 
iSv/6t is nothing but the perturbation added to the respective r.h.s. of the massive Thirring 
model (D). Equations (^)- (|26|) can also be derived via the perturbation theory for the 



Thirring solitons developed in Ref. [50 



From the point of view of the inverse scattering transform method the ansatz (^ contains 
not only the A^-soliton solution but also the contribution of radiation as well. However, due 
to large separations between the solitons in the train, the radiation component is negligible. 
Thus, I can use the adiabatic perturbation theory for derivation of evolution equations for 
the soliton parameters {}a, Va, Xa, and 6a- For the same reason, it is sufficient to consider 
the interaction between the neighboring pulses only (detailed discussion can be found in 
Ref. i). 

First, one should compute the perturbation functions defined in (^) and (PB]). Substi- 
tution of the ansatz (^ into the massive Thirring model (Q) and expansion of the cubic 
terms leads to the following formulae for the perturbation-induced evolutions (I consider the 
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interaction between the neighboring pulses only) 



I— = - {\Ua\'^Vl3 + '^'^(i{UaU*p}Vc 



(27) 



and, consequently, 



.S\Va 
I 



|2 



2i\uo\^\m.{vaVl}, 



(28) 



.5\uc 



/3=QTl 



Before giving the perturbation-induced evolution of the soliton parameters some remarks 
must be maid on the details of the approximation due to the inequalities (0). The r.h.s.'s in 
(^) contain the small parameter e (|12D. I consider the soliton-soliton interaction in the first 
order with respect to e. Hence, due to the presence of the small parameter, the differences 
between the a-th soliton amplitude and rapidity and the average values of these quantities 
(given by (^) are negligible in the terms accounting for the inter-soliton interaction. 

Substitution of (|3) and (|8D into (|2|) and (||) and the result into (|g) and (|l]) gives 
the following equations for the amplitudes and rapidities 

d^9a V- 2sin^(27?) ia , , 

-T^= E ^e-l^-^'sinvl/^^, 29 

^i^^i coshy 

^= E ^^^^^^^4^^eH^-lcosv^.,, (30) 
dt ^i^^, coshy 

where 

A„/3 = sin(2i^) cosh(?/)(x/3 - x^), 

= 5q - - cos(2??) sinh(|/)(a;a - X/j). 

Here exp(— | Aq^jI) = 0{t). It can be easily verified that equations (^91) and ( pO|) do not 
affect the average amplitude d and rapidity y. 

The evolution equations for and 5^ (see (^) and (0)) are comprised of two ad- 
dends, which account for the unperturbed and perturbation-induced evolution (the latter 
contain the soliton overlap parameter e). Hence, one can neglect the terms accounting for 
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the perturbation-induced evolution of these parameters as compared to their unperturbed 
evolution: 

= tanh?/„, = -cos(279„)sech?/„. (31) 

Now everything is ready for derivation of the complex Toda chain. Let us differentiate 
the following quantity — A^/j + i'^ais (/9 = « ± 1): 

^(-A„^ + = sin(27?)sech(y)[(y„ + 2^^) - {yp + 2i^p)l 

where the second order terms in {da — ^) and (y^ — y) are dropped. On the other hand, 
from equations (^) and (|30|) one derives 

X ■s-^ 4sgn(AQ,fl) sin^(2-(9) , .-r m 

—Aya + 2ida)= E I- ^ exp{sgn(A^^)(-A^^ + 2^^^)}, 

at i3=a^i cosii y 

or using the numeration Xa+i — > 0, i.e. sgn(Aa^Q+i) > 0, for the solitons in the train: 

d 4sin^(2'i9) 

— iya + 2ida) = 7—— (exp {-A^^^+i + i'^aa+i} - exp {A^a^i - i'^aa-l}) ■ (32) 

at cosh y 

Introduce an average phase 

6 = — cos(2'(9)sech(|/)t, 
and the following complex variables associated with each soliton 

qa = — sin(2^9)cosh(?/)xa — i[Sa — 5 + avr — cos(2'^) smh.{y)xa] + 2a;ln[2 sin(2'(9)]. 

(The average phase in the formula for qa eliminates the constant phase gradient.) Then 

g±(..±i-<?.) ^ -Asm\2d)exp{TAaa±i ±t^aa±l). (33) 

Differentiating qa and neglecting the second order terms I get 

dqa ^ ^ ^ ^ 

= — sin(2i?) coshy tanh?/a — z[sin(2'(?)sech(|/)(2i?a — 2^9) + cos(2'(9)sech|/tanh(?/)(?/Q, — y) 

— cos(2'i?) sinh y tanh ya] 

(here the average phase 5 plays an important role for the expansion over ■^a—'^)- The second 
differentiation (with removal of the second order terms) gives: 

^ = - sin(2^)sech(y) + 2.^) . (34) 
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At the same time, as it follows from (|32D and ([33|), 

^i-Va - ii^da - 2d)) = sin(2?9)sech(y) (e^^+i"''" - e^^-'^-i) . (35) 
Therefore, what is left is to introduce a new time variable 

r = sin(2'(9)sech(|/)t. 



Then equations ( [5^ ) and ( P^D give the complex Toda chain for the train of Thirring solitons: 

dr^ 

Here it is assumed that Re{go} = and Re{gAr+i} = — oo (i.e., Xq = — oo and xjy+i = oo). 

IV. THE GENERALISED CTC FOR THE TRAIN OF GAP SOLITONS 

In derivation of the complex Toda chain for the optical gap system I will use a one-to-one 
mapping, found recently in Ref. between the optical gap system (|I|) and the following 



generalisation of the massive Thirring model, the 7-system for short, 
t {Vt -V.)+U+ IWI^V + 7~ (1 V|2 - V = 0, 

(36) 

i {Ut +U^) + V+ \V\^U + 7+ - I Vp) W = 0. 



The transformation relating the two systems is as follows. Let U{x,t) and V{x,t) be a 
solution of the 7-system (^6]) with the following 7±: 

" l + pcosh(22/o)' ^^^^ 



then 



I Ei{X,T)\ _ e^^(^'*) |^e-W2v(x,t) \ 

\^E2(X,T) j [l + pcosh(2|/o)]i/2 ^ ey°'^U{x,t) j' 



(38) 



where 



^ X - tanh(yo)r ^ T - tanh(yo)X 
~ [1 - tanh2(|/„)]V2' - [1 _ tanh2(|/„)]V2' 

is a solution to the optical gap system (0) with the phase given by the following system 
of equations (in the light-cone variables 77 = (t + x)/2 and ^ = (t — x)/2) 

|4(o.-og|VP, (40, 
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(Note that the conservation of the number of particles, i. e., 

ensures the compatibihty of the equations for the phase 

Note that the co-ordinates are related via a Lorentz transformation. The mapping 
can be verified by direct substitution into ([1|) via simple calculations with the use of (p7[), 
(H), and iO). 

The presented mapping is valid for arbitrary solutions of the optical gap system. However, 
the importance of the mapping (^7|)-(40) stems from the fact that by choosing the quiescent 
Thirring soliton, 

y — T~/ TTTTi ^ 



cosh(z — I'd) ' cosh(z + id) ' 

z = — sm{2'd)ix — Xo)^ — = — cos(2'i9), 

dt 

which is a solution to the system (^) due to | V| = one can recover the optical gap soliton 
moving with any given velocity V = tanh?/o. In this case one obtains formula (|TB|) for the 
additional phase phase ip. Although the optical gap system is not Lorentz- invariant, still it 
makes sense to call yo "rapidity" of the gap soliton due to the transformation (|37|)-(^0D. 

A train of well separated gap solitons moving with arbitrary average rapidity ?/„ can 
always be represented via the transformation (p7|)-(^) as a train of A^ well separated almost 
quiescent Thirring solitons. Indeed, application of the mapping (^7D-(^) to equation (H) 
with V = V and U = u, under the conditions \ya\ ^ 1 and y = 0, yields the train of A^ 
well separated gap solitons with nearly equal amplitudes and rapidities, where the average 
rapidity is equal to the given arbitrary value yo- Moreover, if one neglects the terms of 
order 0(e), then the additional phase ipa of each gap sohton in the train is determined by 
equation ([T6|) with evident changes: ^ -d, z ^ Za and Xo — ^ Xa- 

Convenience of the 7-system ( |36D for the analytical study of gap solitons is based on the 
two following facts. Firstly, the quiescent Thirring soliton satisfies \U\ = \V\. Hence the last 
terms in the 7-system are small if the solution under study is close to the quiescent Thirring 
soliton, or, in terms of the optical gap system, the solution is close to the gap soliton. 
Secondly, the parameters ■y± are bounded for all values of yo and p (including p = 00). 
Therefore the use of the equivalent 7-system allows one to apply the perturbation theory 
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developed for the train of almost quiescent Thirring solitons to the train of gap solitons for 
arbitrary values of the self-phase modulation parameter p. 

Hence, derivation of the complex Toda chain for the gap soliton train can be done in much 
the same way as the derivation of the complex Toda chain for the train of Thirring solitons 
(more precisely, almost quiescent Thirring solitons). The only difference is that there are 
additional small perturbations given by the terms with 7-1- in (P^D, the 7-terms for short. 
Let us first calculate their contribution to the evolution of the soliton parameters and then 
calculate evolution of qa, defined in a similar was as in the previous section, with account 
of these terms as well. 

Below I will take into account that gap soliton train is transformed by the mapping (|37D - 
(^) with Uo = y to the train of almost quiescent Thirring solitons (in the variables V and 
U). For instance, the latter train has |/ = and \ya\ <^ 1. Consider the a-th soliton in such 
train. From ([36| ) one gets expanding the cubic terms: 



^^ = - E (|Wa|'V^ + 2Re{W„W;}Va)-7-(|V«p-|W„|2)V„ 

- l-[^\Va?Vp + V^^Vl-\U^\^Vp-2Re{UJAl]Vo). 

Let us separate the r.h.s. into two parts. The first part is just the same as in the case of 
the massive Thirring solitons, while the second, given by the following formula 



- 7-(2|Va|V + Va'V;-|W„|\/3-2Re{W,W|}V„), (41) 

/3=QTl 

is due to the 7-perturbation and is specific for the gap soliton train only. 1 will detailly 
consider only the 7-perturbation since the first part is accounted just in the same way as 
for the train of almost quiescent Thirring solitons with the same resulting formulae. 

In formula ( ^T]) the first term on the r.h.s. is due to the self-interaction of the gap soliton 
(due to the 7-terms in (|36|), if the rapidity ?/a 7^ 0), the rest account for the inter-sohton 
interaction in the train. The latter terms contain the small parameter e defined in (p!2[). 
I will use the same approximation which has been used for the derivation of the complex 
Toda chain for the massive Thirring model. Additionally one can neglect the difference 
between the modules of Vq, and Ua when calculating the contribution from the inter-soliton 
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interaction terms. This is because the inter-sohton interaction terms already contain the 
small parameter e and \ya\ <^ 1, thus one can put = there. Then, the "variational" 
derivatives simplify considerably: 



Mir' =-^- 



{W-M)V^- J2 47-sin(^,^)Im{Vo,Wo4V„, (42) 



M^) =-7+(|W„|'-|V„|')v„- 47+sin(^,^)Im{W„,Vo/3}W„, (43) 

where Voa = exp{— i^ajVo,, Uoa = Gxp{—i6a}Ka, and "^a/s = — Sf3- Formula (^) obtains 
from (^) by the evident substitution V ^ U, U ^ V and 7_ —>■ 7+. Note that from 
formulae and (E3D it follows that 



Consider first the contribution to evolution of the soliton parameters coming from the 
inter-soliton interaction 7-terms, i.e., the first terms on the r.h.s.'s in formulae ( ^2]) and (^31). 
First, the perturbation functions given in (p5|) and ( p6D must be calculated. The inter-soliton 
interaction terms give the following contributions to the necessary functions: r\\{za) = and 

(I- \ \- ^. sin^(2^)e-l^"/^lsin(^^^) . 

^r^i [cosh(22;„) + cos(2*')J^ 

X "(7+ - 7_)e^'^ + (7+6^^'^ - 7-e-2^'^)e-^'^] , (44) 

where ka = — exp{— ?/q, — i{}a} and Aq^j = sm{2{}){xf3 — Xa)- (On the r.h.s. of ( ^4|) the 
terms of the second order in -(^q, — -(9 are neglected due to the small multiplier exp(— | Aq/3|) = 
C(e).) Now, substitution of the expression (|^) into (pO[) and (pID leads to the following 
contributions to evolution equations for i^q, and i/a- 



dt 



I =0' 

7 



dt 



I = - ^ 4Ksin^(2i9)e"l'^"''l sin^^^. 



where 

ptanh(2?/) M - sm^M) 



K 



p + sech(2?/) sin^(2^9) 
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What concerns the other two parameters Xa and Sa, similar as in section 3, the inter-sohton 
interaction is of order 0{e) and its contribution to evolution of the soliton position and 
phase can be neglected as compared to their unperturbed evolution. 

Now let us consider the contribution to evolution of the soliton parameters coming from 
the self- interaction 7-terms, i.e., the first terms in (^) and (^5]). I get the following contri- 
butions to the functions in ( P^D and (PE|): 

r\\ (za) = 0, r^{ka, Za) + r^^kl, -Za) = 0, (45) 



[cosh(22;Q,) + cos[2Va)\ '- 

(46) 

(here it is taken into account that terms of the second order in ya are negligible). From ( ^5]) 
it follows that the contributions from the self-interaction 7-terms to evolution equations for 
the amplitude and rapidity i/a vanish, while substitution of (|45|) and (|46|) into (|2^ ) gives 
due to \ya\ ^ 1: 

J„ = 2 {(7_ + 7+)[2 sin2(2^„) - 2^9, sin(4^9„)] + 2(7. - 7+)2^„ sm\2i»^)}y^. 

Thus the contributions from the self-interaction 7-terms to evolution of Xa and 6a are de- 
termined by the following coefficients 

Im{Ja} 
2sin2(2i9,)y, 

I can take just the average values of these coefficients (denoted below as u and /i), because 
the following combinations Uaya and fiaVa will enter the evolution equations for the soliton 
parameters and ^ 1. In other words, one can throw away the terms of the second order 
in — ^ and ya- Taking into account the definition of the 7±, where yo = y, I obtain: 

^ _ 4p(2^9cot(2^9) - 1) _ 4ptanh(2y) - 
p + sech{2y) ' ^ p+sech(2|/) 

Let us collect all the contributions, i.e., the terms same as for the Thirring soliton train 



= -17—277^T\— = (7+ - 7-)2t9a. 



(see equations p9| ) -(|3T| ) ) and those accounting for the 7-terms, and write down the corre- 
sponding evolution equations for the parameters of the a-th gap soliton. They read: 

d^„ 



2sin=^(279)e-l^'^''lsin^„/3, (47) 
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"^^^ ^ 4sin3(2^9)e-l^'^^l[sgn(A,^)cos^«^-«;sm^«;3], (48) 



dt 



d6a 



dt 

dXa 



cos(2^9a) + /i sm(2'^)y«, (49) 
(1 + iy)ya. (50) 



dt 



Here 



= sin(2?9)(a;/3 - x„), ^^/j = 5^ - 5/?. 

It is easy to see that the averages "t} and y are not affected by equations (B^) and 
Equations (^71) -(^) are similar to those for the (almost quiescent) Thirring solitons, however 
there are additional terms in the evolution equations for 6a, and Xa- 

Let us now derive the generalised complex Toda chain corresponding to equations (W^) - 
(^). As the derivation is quite similar to that for the quiescent Thirring solitons I will skip 
some details. As in the case of the massive Thirring model, introduce the modified time 

T = sin(2?9)t, 

an average phase 

5 = - cos{2d)t, 
and the complex variables Qa for each soliton: 

qa = -sm{2{})xa-i{Sa-S + a7f) + 2aln[2sm{2&)]. (51) 

Differentiating and throwing away the second order terms one obtains 

^ = -{(1 + + t^)y^ + t[2^^ - 2^]}. 

dr 

Differentiation of this formula gives 

-r^ = - 4sgn(A„^) sin2(2^9) exp {sgn(A,^)(-A„^ + 2^,^)} 

-{v + i/i)Re 4sgn(A„^) sin^{2§) exp {sgn(A«^)(-A„^ + i^ap)] 
+ z/ + i/i)Im 4sgn(A„^)sin2(279)exp{sgn(Ac.^)(-A„^ + i^«^)}. 
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Taking into account the numeration of the sohtons in the train, which is given by Xa+i —Xa > 
or Aaa+i > 0, and the following identity 



I obtain a generalised complex Toda chain for the train of well separated gap solitons 
with nearly equal amplitudes and rapidities: 

^ = (1 + Ap)(e'?"+i-'?" - e""-'?"-^) + 5^(6"-+!-"" - e^S--?:-!)^ (52) 

where Ap and Bp are p-dependent coefficients: 

Ap = ^{u - K/i + + u) + /i]}, 

Bp = i{z/ + K/i - + z/) - /i]}. 

As usual, Re{go} = 00 and Re{gAr+i} = —00. 

Though in equation (|52|) and in the definition of go, (51) I still have the variables r, Xa, 



Sa, and 6 defined through the co-ordinates x and t (see formula (^)), it is easy to reverse 
to the co-ordinates X and T of the optical gap system (|l|). Indeed, to this end one should 
use the transformation (^) (with i/o = y) for the position Xa and the central phase 5a of 
the gap soliton (the phase at X = X^): 

Xa = cosh.{y)[Xa — tanh(|/)T], (53) 

8^ = [-cos(2?9«) +/isin(2t9)?/Jt = [- cos(2i9c.) + /i sin(2^)7/«]cosh(j/)[T - tanh(^)A:„]. (54) 

Also one must use the time transformation dT = cosh(?/)dt in the definition of r and the 
average phase 5: 

T = sin(2'i9)sech(|/)T, 6 = — cos(2'(9)sech(|/)T. 

Now it is evident that, if equation (|53D is used in the definition of qa (|5lD , the term linear 
in T will not give contribution neither to the difference qa — q/3 nor to the second derivative 
of qa, hence it can be dropped. Further, notice that from the r.h.s. of equation (^) only 
the term linear in T will appear in qa (0), if one simply changes the time: t —>■ sech(|/)T. 
Hence, the term proportional to Xa in (|54D must be subtructed from the central phase 6a- 
In doing so, one can neglect the difference between 'da and d due to the inequalities (^ 
and that evolution of ^a is of order 0{e) (i.e., we throw away the second order terms from 
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the second derivative of qa)- Thus we have arrived precisely at the quantity Qa given by 
equation (0), where the shift of the sohton rapidities is taken into account: Ua y + Ua- 
Therefore the result of section 2 is proven. 



V. COMMENTS 



The complex Toda chain model proves to be an universal model for the adiabatic de- 
scription of the train interaction/propagation of solitons in nonlinear PDEs. Indeed, it was 
shown to describe the train propagation of pulses in the nonlinear PDEs of the whole NLS 
hierarchy (i.e., the PDEs associated with the familiar Zakharov-Shabat spectral problem 
||5T| , |52|). More recently, the complex Toda chain was derived for the soliton train of the 
modified NLS equation [110]. This PDE is associated with the quadratic bundle, also known 



as the Wadati-Konno-Ichikawa spectral problem [^]. 

In this paper, the complex Toda chain is shown to describe the soliton train propagation 
in the massive Thirring model. Note that, as it is mentioned in Ref. ]T0|], the massive 
Thirring model is just another representative of the modified NLS hierarchy. Thus the 
complex Toda chain arises in the adiabatic description of the soliton trains in the hierarchy 
of nonlinear integrable PDEs associated with the quadratic bundle as well. This is in favor 
of the universality of the complex Toda chain. 

In construction of the perturbation theory for the massive Thirring model we have used 
the associated Riemann-Hilbert problem |^9|. The usage of the Riemann-Hilbert problem 
allows one to develop the perturbation theory in an unified way for the entire hierarchy 
(see, for instance, Ref. 0], where this was done for the vector NLS hierarchy). Moreover, 
the perturbation-induced evolution equations for the spectral data have one and the same 
form for all integrable PDEs (one can compare the results of Refs. [^5|, |5^). This gives a 
possibility to prove the universality of the complex Toda chain using the approach based on 
the Riemann-Hilbert problem. This is one of the directions for future work. 



In view of recent experimental observation |Q of the multiple gap soliton formation 
in optical fibers with index grating, it is important to have an analytical approach for 
decription of interaction of optical gap solitons. In this paper an analytical approach is 
developed for the train interaction/propagation of gap solitons: the adiabatic propagation 
of gap solitons with nearly equal amplitudes and velocities is governed by a generalised 
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complex Toda chain with nodes. 

Here I should point out that, due to non-integrability of the optical gap system, the train 
of gap solitons may become unstable. This instability is the result of the soliton-radiation 
interaction and is beyond the adiabatic approximation. However, the gap soliton is stable 
if the soliton amplitude lies below the instability threshold (see for details Ref. In 
that case, the generalised complex Toda chain ( [T5| ) can be applied. Though in accordance 
with the Ref. the generalised complex Toda chain is not integrable, it is just a finite di- 
mensional dynamical system and can be investigated by the standard techniques. Moreover, 
in accordance with discussion of Ref. , one can systematically include various additional 
perturbations of the optical gap system into the complex Toda chain. This is the direction 
for further work. 
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